The setup of the tile part imitated the assignment I did for STA303.

Codes are adapted from BCB420 lecture slides.

The journal page for A3 is: https://github.com/bcb420-2022/Xinyi_Xu/wiki/12.-Data-set-Pathway-and-Network-Analysis

1 Introduction

We obtained the raw data from article: Transcriptomic Similarities and Differences in Host Response between SARS-CoV-2 and Other Viral Infection, with GEO(Barrett et al. 2012) id GSE152641(Thair et al. 2021). The data contains 20460 observations with 87 variables, there are in total 86 participants, among them, 62 are COVID-19 patients and 24 are healthy individuals.

We created our normalized data set by removing duplicated and unwanted observations. Box plots and density plots were used to visualize the data and then we perform TMM normalization on our primary cleaned data, after doing so, we noticed that there is a significant difference between normalized and un-normalized data. The entrez_id used by raw data set was transferred into the hgnc_symbol or corresponding gene name for the normalized data set and duplicates were removed again. Our normalized data set has 14413 observations.

Based on the normalized data, we used package limma(Ritchie et al. 2015) and edgeR(Robinson, McCarthy, and Smyth 2010) (McCarthy, Chen, and Smyth 2012) to propose models and calculate the corresponding p-values. Volcano plots and heat maps are graphed for each model. We used a threshold 0.05 to identify genes that are significant, and filter out the unwanted gene names. Three txt files containing upregulate, downregulate, and ranked genes respectively were created for further analysis.

Note the setup section which load packages are hided.

# We load the ranked data we had created in Assignment 2 and skim it. 
ranked_data <- read.table(file="covid_rankedgenelist.rnk",
                                    header = TRUE,sep = "\t",
                                    stringsAsFactors = FALSE,
                                    check.names=FALSE)
kable(ranked_data[1:5,1:2], type="html", caption = "Table1. Skim on the ranked data")
Table1. Skim on the ranked data
Genename rank
H2BC6 -16.14800
POLQ -16.14481
APOBEC3A -15.45742
APOBEC3A_B -15.45742
CC2D2A -14.80352

The data set that we use for GSEA is the ranked data that we have created in Assignment 2, this data set includes all the ranked values and no threshold was applied, thus it has 14413 observations and 2 variables, the number of observations is the same as our normalized data set we have used in Assignment 2, detailed codes for creation of ranked data is in Thresholded over-representation analysis section of Assignment 2 and was slightly modified and documented again on journal page for A3.

2 Non-thresholded Gene set Enrichment Analysis

2.1 What method did you use? What genesets did you use? Make sure to specify versions and cite your methods.

We used GSEA pre-ranked analysis. The specific values and versions we used are:

Version: 4.2.3

Enrichment statistic: weighted

Max size: 200

Min size: 15

Collapse: No_Collapse

The genesets we used was adapted from the baderlab geneset(Merico et al. 2010a) collection from March 1, 2021 with GOBP all pathways and no GO IEA. (Codes are documented in journal GSEA)

2.2 Summarize your enrichment results.

Note the idea of adding tables come from: https://bookdown.org/yihui/rmarkdown-cookbook/kable.html

                                       Table2. Layout of the Enrichment result
Stats na_pos na_neg
Upregulated 2205 3088
FDR < 25% 323 1092
significantly enriched at nominal pvalue < 1% 189 558
significantly enriched at nominal pvalue < 5% 376 910
Top Gene sets EUKARYOTIC TRANSLATION INITIATION%REACTOME DATABASE ID RELEASE 79%72613 HALLMARK_INTERFERON_ALPHA_RESPONSE%MSIGDB_C2%HALLMARK_INTERFERON_ALPHA_RESPONSE

In total 5293 gene sets are used. na_pos gene sets are genes that are upregulated in healthy human individuals, and na_neg gene sets are genes that are downregulated in healthy human individuals.

                                       Table3. Top Genes Sets for na_pos and na_neg
Stats EUKARYOTIC TRANSLATION INITIATION%REACTOME DATABASE ID RELEASE 79%72613 HALLMARK_INTERFERON_ALPHA_RESPONSE%MSIGDB_C2%HALLMARK_INTERFERON_ALPHA_RESPONSE
Upregulated in class na_pos na_neg
Enrichment Score (ES) 0.6688145 -0.8105112
Normalized Enrichment Score (NES) 2.506432 -2.9652746
Nominal p-value 0.0 0.0
FDR q-value 0.0 0.0
FWER p-Value 0.0 0.0
Leading-edge 81 58
Top Gene EIF4B EPSTI1

2.3 How do these results compare to the results from the thresholded analysis in Assignment #2. Compare qualitatively. Is this a straight forward comparison? Why or why not?

It is not straight forward comparison, because the number of genes sets involved in GSEA and thresholded analysis is different, the total gene sets involved in thresholded analysis is 7802, while after filtering, the total GSEA gene sets used are 5293.

But the results we get from thresholded analysis and GSEA are very similar. Both of the top two gene sets we get from GSEA are among the top gene sets we get from REACTOME in thresholded analysis. And by observing the top gene sets of both thresholded analysis and GSEA, they do overlap with each other, like HALLMARK_INTERFERON_GAMMA_RESPONSE%MSIGDB_C2%HALLMARK_INTERFERON_GAMMA_RESPONSE, and GTP HYDROLYSIS AND JOINING OF THE 60S RIBOSOMAL SUBUNIT%REACTOME%R-HSA-72706.2 which are the second top gene sets given by GSEA for na_neg and na_pos, are also among the top gene sets of thresholded analysis, similar examples would be: DEFENSE RESPONSE TO BACTERIUM%GOBP%GO:0042742, CYTOPLASMIC TRANSLATION%GOBP%GO:0002181 and so on.

3 Visualize your Gene set Enrichment Analysis in Cytoscape

We use Cytoscape(Shannon et al. 2003) version 3.9.1 for MaxOS. Apps like EnrichmentMap(Merico et al. 2010b), and AutoAnnotate(Kucera et al. 2016) are installed from Cytoscape app store.

3.1 How many nodes and how many edges in the resulting map? What thresholds were used to create this map? Make sure to record all thresholds. Include a screenshot of your network prior to manual layout.

There are in total 534 nodes and 2771 edges. The thresholds used is listed in the figure below: The FDR q-value cutoff was set to 0.05, and all other values remained default.

Fig1. Thresholds used to create the Enrichment Map

Fig2. Network prior to manual layout

Fig3. Legend

3.2 What parameters did you use to annotate the network. If you are using the default parameters make sure to list them as well.

The plugin package used to annotate the network is: AutoAnnotate(Kucera et al. 2016). We used the default choice plus the Layout network to prevent cluster overlap.

Fig4. Annotated Network without prevent cluster overlap

Fig5. Annotated Network without prevent cluster overlap

Fig6. Annotated Network with prevent cluster overlap

Fig7. Legend

3.3 Figure

We use the publication ready button in the Enrichment map section.

Fig8. Publication ready figure

3.4 Collapse your network to a theme network. What are the major themes present in this analysis? Do they fit with the model? Are there any novel pathways or themes?

Fig9. The Collapsed Network with no overlap

Fig10. The Collapsed Network with overlap

There are in total 173 clusters and 41 edges in both Networks. The major themes present in this analysis are: cell cycle, translation and transcription, immune response, viral life cycle, cancer, and DNA repair.

We believe that the result we get matches with the model. As EUKARYOTIC TRANSLATION INITIATION%REACTOME DATABASE ID RELEASE 79%72613 is clearly gene set that is involved in translation and transcription, while HALLMARK_INTERFERON_ALPHA_RESPONSE%MSIGDB_C2%HALLMARK_INTERFERON_ALPHA_RESPONSE is a gene set that involved in immune system, and interferons are proteins which indicate that there are invaders or cancer cells in the body, and this matches with our purpose of investigation as we study the difference between healthy and Covid19 individuals. This result of GSEA is matched with our theme network, there is no clear novel pathways or themes.

4 Interpretation and detailed view of results

4.1 Do the enrichment results support conclusions or mechanism discussed in the original paper? How do these results differ from the results you got from Assignment #2 thresholded methods?

The enrichment results support the article: Transcriptomic Similarities and Differences in Host Response between SARS-CoV-2 and Other Viral Infection(Thair et al. 2021), as the article focused on covid 19 and other viral infections which matches with our theme: viral life cycle and immune response. Immune response and cancer was also discussed in the article where the authors claim that neutrophil-to-lymphocyte ratios are used to mark cancer severity(Diao et al. 2020) (Liu et al. 2020) (Lagunas-Rangel 2020) (Qin et al. 2020). The article also stated that translation pathway like translation of ferritin is affected with COVID19 infections. Moreover, the article mentioned that genes involved in neutrophil activation, innate immune response, immune response to viral infection, type-I interferon signaling, cytokine production and lymphocyte differentiation and T-cell activation and regulation(Thair et al. 2021) differ between healthy and Covid19 individuals, which highly matches with our enrichment results.

Novel pathways is DNA repair (notch expression ercc6) which was not discussed in the article.

The overall results matches with our Assignment #2 thresholded analysis results, as in our analysis in Assignment 2, we found most genes involved in the pathways mentioned above are among the top ranked genes in g:Profiler(Reimand et al. 2007) analysis. The results are consistent and do not show significant difference. The slight difference may due to the fact that the total gene sets involved in thresholded analysis is 7802, while after filtering, the total GSEA gene sets used are 5293 as we have discussed in the previous section.

4.2 Can you find evidence, i.e. publications, to support some of the results that you see. How does this evidence support your result?

Evidence and discussions that we did in Assignment #2 also supports our enrichment results. Besides that, several articles(Diao et al. 2020) (Liu et al. 2020) (Lagunas-Rangel 2020) (Qin et al. 2020) mentioned in the previous section, provided an evidence to supporting our result: immune response and cancer are among the major themes present in this analysis.

The article(Pánico, Ostrosky-Wegman, and Salazar 2022) claims that viral infections could lead to DNA damage and might also increase the risk to cancer development, which supports our finding that there are major themes: DNA repair and cancer involved in the major themes, though DNA repair is not mentioned in our original paper.

4.3 Add a post analysis to your main network using specific transcription factors, microRNAs or drugs. Include the reason why you chose the specific miRs, TFs or drugs (i.e publications indicating that they might be related to your model). What does this post analysis show?

I choose to do post analysis on drugs.

# Codes imitated the GSEA section
# We first download the drugbank
gmt_url = "https://download.baderlab.org/EM_Genesets/current_release/Human/symbol/DrugTargets/"
# list all the files on the server
filenames = getURL(gmt_url)
tc = textConnection(filenames)
contents = readLines(tc)
close(tc)
# get the gmt that has all the pathways and does not include terms inferred
# Human_DrugBank_approved_symbol.gmt  
rx = gregexpr("(?<=<a href=\")(.*.approved.*.)(.gmt)(?=\">)", contents,
    perl = TRUE)
gmt_file = unlist(regmatches(contents, rx))
dest_gmt_file <- file.path(gmt_file)
download.file(paste(gmt_url, gmt_file, sep = ""), destfile = dest_gmt_file)

4.3.1 Reason

In the original paper(Thair et al. 2021), the author mentioned that drug pegylated interferon alpha was approved by FDA using to treat HCV virus infection and previously showed a promising result when combining with ribavirin to treat MERS. In the references used by the original paper, we noticed that most references used by (Tan et al. 2004) (The reference of the original paper) was related to interferon alpha2b and some of them even considered the combination of ribavirin, while (Nile et al. 2020) mentioned that for MERS, usage of oral ribavirin combined with subcutaneous pegylated IFN alpha-2a resulted in noticeable improvement in the chance of survival. Thus, it would be a good idea for us to take both pegylated IFN alpha-2a and pegylated IFN alpha-2b together with ribavirin in to account while doing our post analysis, as according to our original paper, the drug has mechanisms that may be helpful in treating COVID19 (also mentioned in (Nile et al. 2020)).

Note, we used Mann-Whitney(one sided less) Test with a cutoff at 0.05.

4.3.2 Known Signitures

The following is the value of the three drugs we get:

Fig11. The Known Signitures

However, after analysis, we very unfortunately found no edges that pass the signature cutoff, and even when we expand the network, still no edges was found. This may suggests that pegylated IFN alpha-2a and pegylated IFN alpha-2b together with ribavirin may not show a similar promising result when being used to treat COVID19 as it did while being used to treat MERS.

Thus we would use Exploratory post analysis to see if any other drugs has passed the cutoff.

4.3.3 Exploratory

Fig12. Top Drugs that passed the Cutoff

We noticed that FOSTAMATINIB, GLUTAMIC ACID and PRALSETINIB are among the top drugs, which indicates that they might be the possible drug candidates being used to treat COVID19. As we explore in Google Scholar(Vine 2006), we noticed that Strich(Strich, Ramos-Benitez, et al. 2021) (Strich, Tian, et al. 2021) had done a randomized trial on FOSTAMATINIB and claimed that FOSTAMATINIB shows a promising result in saving individuals who has severe disease symptoms.

While Naringenin, a drug that interact with GLUTAMIC ACID, was investigated in 2020 to treat COVID19(Tutunchi et al. 2020), and PRALSETINIB was approved by FDA to treat lung and thyroid cancers(Kim et al. 2021). However, no published paper clearly indicated any promising results from clinical trials of these two drugs when being used to treat COVID19.

And again, when we use the theme network, very unfortunately, we found no edges that pass the signature cutoff. Indicates that the drugs with the approved symbol from DrugBank(Wishart et al. 2018) may not yield a significant result while being used to treat COVID19.

However, when we expand all, we are able to witness a variety of edges linked between the nodes of our Enrichment map and the three top drugs. Among them, FOSTAMATINIB exhibits the most linkages, which indicates that the drug is a possible candidate for treating COVID19.

Fig13. Top Drugs With Expanded Enrichment Map (whole graph)

Fig14. Top Drugs With Expanded Enrichment Map (close skim)

References

Bankar, Renuka, Kruthi Suvarna, Saicharan Ghantasala, Arghya Banerjee, Deeptarup Biswas, Manisha Choudhury, Viswanthram Palanivel, et al. 2021. “Proteomic Investigation Reveals Dominant Alterations of Neutrophil Degranulation and mRNA Translation Pathways in Patients with COVID-19.” iScience 24 (3): 102135. https://doi.org/https://doi.org/10.1016/j.isci.2021.102135.
Barrett, Tanya, Stephen E. Wilhite, Pierre Ledoux, Carlos Evangelista, Irene F. Kim, Maxim Tomashevsky, Kimberly A. Marshall, et al. 2012. NCBI GEO: archive for functional genomics data sets—update.” Nucleic Acids Research 41 (D1): D991–95. https://doi.org/10.1093/nar/gks1193.
Blighe, Kevin, Sharmila Rana, and Myles Lewis. 2020. EnhancedVolcano: Publication-Ready Volcano Plots with Enhanced Colouring Andlabeling. https://github.com/kevinblighe/EnhancedVolcano.
Davis, Sean, and Paul Meltzer. 2007. “GEOquery: A Bridge Between the Gene Expression Omnibus (GEO) and BioConductor.” Bioinformatics 14: 1846–47.
Diao, Bo, Chenhui Wang, Yingjun Tan, Xiewan Chen, Ying Liu, Lifen Ning, Li Chen, et al. 2020. “Reduction and Functional Exhaustion of t Cells in Patients with Coronavirus Disease 2019 (COVID-19).” Frontiers in Immunology 11. https://doi.org/10.3389/fimmu.2020.00827.
Durinck, Steffen, Yves Moreau, Arek Kasprzyk, Sean Davis, Bart De Moor, Alvis Brazma, and Wolfgang Huber. 2005. “BioMart and Bioconductor: A Powerful Link Between Biological Databases and Microarray Data Analysis.” Bioinformatics 21: 3439–40.
Gentleman, R. 2020. Annotate: Annotation for Microarrays.
Gu, Zuguang, Roland Eils, and Matthias Schlesner. 2016. “Complex Heatmaps Reveal Patterns and Correlations in Multidimensional Genomic Data.” Bioinformatics.
Gu, Zuguang, Lei Gu, Roland Eils, Matthias Schlesner, and Benedikt Brors. 2014. “Circlize Implements and Enhances Circular Visualization in r.” Bioinformatics 30: 2811–12.
Kim, Janice, Diana Bradford, Erin Larkins, Lee H Pai-Scherf, Somak Chatterjee, Pallavi S Mishra-Kalyani, Emily Wearne, et al. 2021. “FDA Approval Summary: Pralsetinib for the Treatment of Lung and Thyroid Cancers with RET Gene Mutations or Fusions.” Clinical Cancer Research 27 (20): 5452–56.
Kolberg, Liis, Uku Raudvere, Ivan Kuzmin, Jaak Vilo, and Hedi Peterson. 2020. “Gprofiler2– an r Package for Gene List Functional Enrichment Analysis and Namespace Conversion Toolset g:profiler.” F1000Research 9 (ELIXIR) (709).
Kucera, Mike, Ruth Isserlin, Arkady Arkhangorodsky, and Gary D Bader. 2016. “AutoAnnotate: A Cytoscape App for Summarizing Networks with Semantic Annotations.” F1000Research 5.
Lagunas-Rangel, Francisco Alejandro. 2020. “Neutrophil-to-Lymphocyte Ratio and Lymphocyte-to-c-Reactive Protein Ratio in Patients with Severe Coronavirus Disease 2019 (COVID-19): A Meta-Analysis.” Journal of Medical Virology.
Liu, Jingyuan, Yao Liu, Pan Xiang, Lin Pu, Haofeng Xiong, Chuansheng Li, Ming Zhang, et al. 2020. “Neutrophil-to-Lymphocyte Ratio Predicts Critical Illness Patients with 2019 Coronavirus Disease in the Early Stage.” Journal of Translational Medicine 18 (1): 1–12.
Long, Quan-Xin, Xiao-Jun Tang, Qiu-Lin Shi, Qin Li, Hai-Jun Deng, Jun Yuan, Jie-Li Hu, et al. 2020. “Clinical and Immunological Assessment of Asymptomatic SARS-CoV-2 Infections.” Nature Medicine 26 (8): 1200–1204.
McCarthy, Davis J, Yunshun Chen, and Gordon K Smyth. 2012. “Differential Expression Analysis of Multifactor RNA-Seq Experiments with Respect to Biological Variation.” Nucleic Acids Research 40 (10): 4288–97. https://doi.org/10.1093/nar/gks042.
Merico, Daniele, Ruth Isserlin, Oliver Stueker, Andrew Emili, and Gary D Bader. 2010a. “Enrichment Map: A Network-Based Method for Gene-Set Enrichment Visualization and Interpretation.” PloS One 5 (11): e13984.
———. 2010b. “Enrichment Map: A Network-Based Method for Gene-Set Enrichment Visualization and Interpretation.” PloS One 5 (11): e13984.
Montojo, Jason, Khalid Zuberi, Harold Rodriguez, Farzana Kazi, George Wright, Sylva L Donaldson, Quaid Morris, and Gary D Bader. 2010. “GeneMANIA Cytoscape Plugin: Fast Gene Function Predictions on the Desktop.” Bioinformatics 26 (22): 2927–28.
Mootha, Vamsi K, Cecilia M Lindgren, Karl-Fredrik Eriksson, Aravind Subramanian, Smita Sihag, Joseph Lehar, Pere Puigserver, et al. 2003. “PGC-1\(\alpha\)-Responsive Genes Involved in Oxidative Phosphorylation Are Coordinately Downregulated in Human Diabetes.” Nature Genetics 34 (3): 267–73.
Morgan, Martin. 2021. BiocManager: Access the Bioconductor Project Package Repository. https://CRAN.R-project.org/package=BiocManager.
Nile, Shivraj Hariram, Arti Nile, Jiayin Qiu, Lin Li, Xu Jia, and Guoyin Kai. 2020. “COVID-19: Pathogenesis, Cytokine Storm and Therapeutic Potential of Interferons.” Cytokine & Growth Factor Reviews 53: 66–70.
Pagès, Hervé, Marc Carlson, Seth Falcon, and Nianhua Li. 2020. AnnotationDbi: Manipulation of SQLite-Based Annotations in Bioconductor. https://bioconductor.org/packages/AnnotationDbi.
Pánico, Pablo, Patricia Ostrosky-Wegman, and Ana María Salazar. 2022. “The Potential Role of COVID-19 in the Induction of DNA Damage.” Mutation Research/Reviews in Mutation Research 789: 108411. https://doi.org/https://doi.org/10.1016/j.mrrev.2022.108411.
Qin, Chuan, Luoqi Zhou, Ziwei Hu, Shuoqi Zhang, Sheng Yang, Yu Tao, Cuihong Xie, et al. 2020. “Dysregulation of Immune Response in Patients with Coronavirus 2019 (COVID-19) in Wuhan, China.” Clinical Infectious Diseases 71 (15): 762–68.
R Core Team. 2020. R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing. https://www.R-project.org/.
Reimand, Jüri, Meelis Kull, Hedi Peterson, Jaanus Hansen, and Jaak Vilo. 2007. g:Profiler—a web-based toolset for functional profiling of gene lists from large-scale experiments.” Nucleic Acids Research 35 (suppl_2): W193–200. https://doi.org/10.1093/nar/gkm226.
Ritchie, Matthew E, Belinda Phipson, Di Wu, Yifang Hu, Charity W Law, Wei Shi, and Gordon K Smyth. 2015. limma Powers Differential Expression Analyses for RNA-Sequencing and Microarray Studies.” Nucleic Acids Research 43 (7): e47. https://doi.org/10.1093/nar/gkv007.
Robinson, Mark D, Davis J McCarthy, and Gordon K Smyth. 2010. “edgeR: A Bioconductor Package for Differential Expression Analysis of Digital Gene Expression Data.” Bioinformatics 26 (1): 139–40. https://doi.org/10.1093/bioinformatics/btp616.
Schulte-Schrepping, Jonas, Nico Reusch, Daniela Paclik, Kevin Baßler, Stephan Schlickeiser, Bowen Zhang, Benjamin Krämer, et al. 2020. “Severe COVID-19 Is Marked by a Dysregulated Myeloid Cell Compartment.” Cell 182 (6): 1419–1440.e23. https://doi.org/https://doi.org/10.1016/j.cell.2020.08.001.
Shannon, Paul, Andrew Markiel, Owen Ozier, Nitin S Baliga, Jonathan T Wang, Daniel Ramage, Nada Amin, Benno Schwikowski, and Trey Ideker. 2003. “Cytoscape: A Software Environment for Integrated Models of Biomolecular Interaction Networks.” Genome Research 13 (11): 2498–2504.
Stempor, Przemyslaw, and Julie Ahringer. 2014. “SeqPlots - Interactive Software for Exploratory Data Analyses, Pattern Discovery and Visualization in Genomics [Version 1; Referees: 2 Approved, 1 Approved with Reservations].” Wellcome Open Research 1 (14). https://doi.org/10.12688/wellcomeopenres.10004.1.
Strich, Jeffrey R, Marcos J Ramos-Benitez, Davide Randazzo, Sydney R Stein, Ashley Babyak, Richard T Davey, Anthony F Suffredini, Richard W Childs, and Daniel S Chertow. 2021. “Fostamatinib Inhibits Neutrophils Extracellular Traps Induced by COVID-19 Patient Plasma: A Potential Therapeutic.” The Journal of Infectious Diseases 223 (6): 981–84.
Strich, Jeffrey R, Xin Tian, Mohamed Samour, Christopher S King, Oksana Shlobin, Robert Reger, Jonathan Cohen, et al. 2021. “Fostamatinib for the Treatment of Hospitalized Adults with COVD-19 a Randomized Trial.” Clin. Infect. Dis.
Subramanian, Aravind, Pablo Tamayo, Vamsi K. Mootha, Sayan Mukherjee, Benjamin L. Ebert, Michael A. Gillette, Amanda Paulovich, et al. 2005. “Gene Set Enrichment Analysis: A Knowledge-Based Approach for Interpreting Genome-Wide Expression Profiles.” Proceedings of the National Academy of Sciences 102 (43): 15545–50. https://doi.org/10.1073/pnas.0506580102.
Tan, Emily LC, Eng Eong Ooi, Chin-Yo Lin, Hwee Cheng Tan, Ai Ee Ling, Bing Lim, and Lawrence W Stanton. 2004. “Inhibition of SARS Coronavirus Infection in Vitro with Clinically Approved Antiviral Drugs.” Emerging Infectious Diseases 10 (4): 581.
Temple Lang, Duncan. 2022. RCurl: General Network (HTTP/FTP/...) Client Interface for r. https://CRAN.R-project.org/package=RCurl.
Thair, Simone A., Yudong D. He, Yehudit Hasin-Brumshtein, Suraj Sakaram, Rushika Pandya, Jiaying Toh, David Rawling, et al. 2021. “Transcriptomic Similarities and Differences in Host Response Between SARS-CoV-2 and Other Viral Infections.” iScience 24 (1): 101947. https://doi.org/https://doi.org/10.1016/j.isci.2020.101947.
Tutunchi, Helda, Fatemeh Naeini, Alireza Ostadrahimi, and Mohammad Javad Hosseinzadeh-Attar. 2020. “Naringenin, a Flavanone with Antiviral and Anti-Inflammatory Effects: A Promising Treatment Strategy Against COVID-19.” Phytotherapy Research 34 (12): 3137–47.
Vine, Rita. 2006. “Google Scholar.” Journal of the Medical Library Association 94 (1): 97.
Wickham, Hadley. 2016. Ggplot2: Elegant Graphics for Data Analysis. Springer-Verlag New York. https://ggplot2.tidyverse.org.
Wickham, Hadley, Mara Averick, Jennifer Bryan, Winston Chang, Lucy D’Agostino McGowan, Romain François, Garrett Grolemund, et al. 2019. “Welcome to the tidyverse.” Journal of Open Source Software 4 (43): 1686. https://doi.org/10.21105/joss.01686.
Wishart, David S, Yannick D Feunang, An C Guo, Elvis J Lo, Ana Marcu, Jason R Grant, Tanvir Sajed, et al. 2018. “DrugBank 5.0: A Major Update to the DrugBank Database for 2018.” Nucleic Acids Research 46 (D1): D1074–82.
Wu, Guanming, Xin Feng, and Lincoln Stein. 2010. “A Human Functional Protein Interaction Network and Its Application to Cancer Data Analysis.” Genome Biology 11 (5): 1–23.
Xie, Yihui. 2014. “Knitr: A Comprehensive Tool for Reproducible Research in R.” In Implementing Reproducible Computational Research, edited by Victoria Stodden, Friedrich Leisch, and Roger D. Peng. Chapman; Hall/CRC. http://www.crcpress.com/product/isbn/9781466561595.